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ABSTRACT 

We report the discovery of an X-ray pulsar in the young, massive Galactic star cluster Westerlund 1. 
We detected a coherent signal from the brightest X-ray source in the cluster, CXO J164710. 2-455216, 
during two Chandra observations on 2005 May 22 and June 18. The period of the pulsar is 10.6107(1) 
s. We place an upper limit to the period derivative of P < 2 x 10~ 10 s s" 1 , which implies that 
the spin-down luminosity is E < 3 x 10 33 erg s _1 . The X-ray luminosity of the pulsar is Lx ~ 
3^2° x 10 33 (D/5 kpc) 2 erg s~\ and the spectrum can be described by a kT — 0.61 tp'^ keV blackbody 
with a radius of i?bb = 0.27 ± 0.03(£>/5 kpc) km. Deep infrared observations reveal no counterpart 
with K<18.5, ruling out a binary companion with M>1M©. Taken together, the properties of 
the pulsar indicate that it is a magnetar. The rarity of slow X-ray pulsars and the position of 
CXO J164710. 2-455216 only 1.6' from the core of Westerlund 1 indicates that it is a member of the 
cluster with >99.97% confidence. Westerlund 1 contains 07V stars with initial masses Mi«35M and 
>50 post-main-sequence stars that indicate the cluster is 4±1 Myr old. Therefore, the progenitor 
to this pulsar had an initial mass Mi>40M Q . This is the most secure result among a handful of 
observational limits to the masses of the progenitors to neutron stars. 

Subject headings: X-rays: stars — neutron stars — open clusters and associations: individual (West- 
erlund 1) 



1. INTRODUCTION 

Most of our knowledge about the masses of the pro- 
genitors to ne utron stars is bas ed on theoretical cal- 
culations fe.g.. iHeger et al1l2003lK Quantitative obser- 
vational constraints are difficult to establish. The few 
previous, tentative estimates have relied on inferring 
these masses from traces of the int eractions of the pro- 
genitors with their surroundings llNqmoto et alJ Il982t 
iMacAlpine et all 119891 iGaensler et alJl2005|) . o~ devel- 
oping scenarios by which indi vidual accreting X-ray bi- 
narie s could have formed (e.g..lErgma fc van den Heuvell 
119981 iWellstein fc Langeil 11999(1 . or on demonstrating 
that the progenitor was a member of a popu lation of co- 
eval stars with well-determined masses (e.g.. lFuchs et alJ 
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1999; Vrh a et al .1 12000: Ei ger et a,lJl2005l IPellizza et all 

2005). The first two methods remain uncertain be- 
cause of their dependence on assumptions in the rel- 
evant models. The third class of results are poten- 
tially the most reliable, although there is sometimes de- 
bate over the whether the neutr on stars and the stel- 
lar associations are related (e.g.. [ Cameron et alJ 120051: 
iMcClu rc-Griffiths fc Gaenslerll2005lK 

Here we report the discovery of a neutron star in Chan- 
dra observations of the Galactic star cluster Westerlund 
1. The cluster contains an exceptional population of 
>50 massi ve, post-main-sequence stars that ar e only 4±1 
Myr old ( Westcrlundl ll987l: iClark et al.l l2*005) . Its total 
mass of >10 5 M Q is contained in a region only «9 pc 
across, which suggests that the stars were born i n a sin- 
gle episode of star formation l)Clark et al J 120 05). This 
makes it an ideal cluster for placing limits on the initial 
masses of progenitors to compact objects. 

2. OBSERVATIONS AND DATA ANALYSIS 

We observed Westerlund 1 with the Chandra X-ray 
Observatory A dvanced CCD Spectrom eter Spectroscopic 
array fACIS-S lWeisskopf et all {2002) on two occasions: 
2005 May 22 for 18 ks (sequence 5411) and 2005 June 18 
for 42 ks (sequence 6283). The event lists were reduced 
using the standard tools and techniques described on the 
web site of the Chandra X-ray Center (CXC). 

We used the wav elet-based algorithm wavdetect 
ijFreeman et all 120021) to identify 238 individual point- 
like X-ray sources in the entire image (which covered 
approximately four times the area displayed in Fig. |JJ. 
Half of these sources are located in the central por- 
tion of the image shown. Above our completeness limit 
of 6 x 10~ 7 ph cm" 2 s" 1 (0.5-8.0 keV), we find that 
the central surface density of X-ray sources is 33l}g 
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Fig. 1. — The central 5' x5' of the image we obtained of Wester- 
lund 1 with Chandra. The image is a composite of two exposures, 
taken on 2005 May 22 and June 18. The count image has been cor- 
rected to account for the varying exposure across the image. We 
have indicated the centroid and core-radius of the cluster derived 
from the locations of the X-ray sources using the dashed line. The 
location of the pulsar is indicated with the solid circle. 
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Fig. 2. — The binned radial distribution of point sources brighter 
than our 50% completeness limit. The distribution is computing 
assuming that its centroid is at a,S = 16 47 03.7, -45 51 00 (J2000). 
The solid line is the best-fit Lorentzian model of the unbinned 
distribution, which has a width of 8o = 0.4 ±0.1 arcmin and a 
central surface density of po = 33^ 



arcmin 2 . For comparison, observations of the Galac- 
tic plane at I = 28° and b = 0.2° reveal only 0.17±0.04 
sourc es arcmin -2 at this flux level (e.g., lEbisawa et alJ 
1200 The surface density of sources can be modeled 
as a Lorentzian function with half- width 0.4±0.1 arcmin 
(Fig. . Therefore, even without considering the proper- 
ties of individual X-ray sources, any one located within 
1'.7 of the core of Westerlund 1 has a >90% chance of 
being a cluster member. 

Most of the X-ray sources have Wolf-Rayet or O 
(WR/O) stars as counterparts in optical or infrared im- 
ages, and are probably colliding- wind binaries. However, 
some of the WR/O stars could be in binaries with accret- 
ing compact objects, and X-ray sources without obvious 
stellar counterparts could be isolated pulsars. We did 
not find any obvious candidates for accreting black holes 
(J. S. Clark et al., in prep). 

The best means of identifying neutron stars is to search 
for rotationally-modulated X-ray emission. Therefore, 
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Fig. 3.— The periodogram of CXO J164710.2-455216, which 
is the only source to exhibit a significant periodic signal. The 
dashed lines denote the powers above which there is a <1% chance 
that noise would produce signals that large. The signal from 
CXO J164710.2-455216 produced Z\ = 58.3 from 398 counts in 
the May observation, and Zj = 135 from 857 events in the June 
observation. The periodicity was significant at the >8<r level. The 
profile of the signal was sinusoidal, and the fractional rms ampli- 
tudes were 53.0(1)% and 54.91(3)%. 



for the brightest sources, we adjusted the arrival times of 
their photons to the Solar System barycenter and com- 
puted Fo urier periodograms using the Rayleigh statis- 
tic (Z 2 ; iBucceri et aTT 119831) . The individual X-ray 
events were recorded with a time resolution of 3.2 s, 
so the Nyquist frequency was «0.15 Hz, which repre- 
sents the limit below which our sensitivity could be well- 
characterised. However, we computed the periodogram 
using a maximum frequency of ss0.6 Hz, to take advan- 
tage of the limited sensitivity to higher frequency signals, 
and to ensure that any observed signal was not an alias. 
We searched the two observations separately, so that a 
search for signals from a single source required 3221 in- 
dependent trials for the May observation, and 6239 tri- 
als for the June observation. The corresponding single- 
source detection threshold powers for 99% confidence 
in each observation are Z\=Yl.l and 13.3, respectively, 
where the powers have been normalised so that Poisson 
noise produces power with a mean value of 1. This power 
can be related to the root-mean-squared (rms) ampli- 
tude of a sinusoidal signal by A — (2Z 2 /iV 7 ) 1 / 2 , where 
iV 7 is the number of photons from the source. A fully- 
modulated signal would have an amplitude of 0.71, so the 
minimum number of counts required to detect a signal 
are 7V 7 =51 and 53, respectively. We searched 8 sources 
above this count limit in the May observation, and 16 
sources in the June observation. 

One periodic signal from the brightest X-ray source 
in the field, CXO J164710.2-455216 (a, 6 = 251.79250,- 
45.87136 [J2000], ±(y.'3 with 90% conf.), significantly ex- 
ceeded the expected noise in both observations (Fig.|2J). 
We refined an initial estimate of the period by com- 
puting pulse profiles from non-overlapping 5000 s in- 
tervals during each observation, and modelling the dif- 
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Fig. 4.— The X-ray spectrum of CXO J164710.2-455216. The 
top panel displays the spectrum in detector counts as a function 
of energy, so the shape of the source spectrum is convolved with 
the detector response. Several models reproduced the data equally 
well (see text). The bottom panel contains the difference between 
the data and the best-fit black body model, divided by the Poisson 
uncertainty on the data points. 
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Fig. 5. — Infrared image a 10"xl0" field around 
CXO J164710.2-455216. The image was taken with the ESO/NTT 
in the K 3 band, and is composed of 10 X 1.2s images obtained on 
2003 June 19 (P.I.: J. Alves). The circle denotes the 0'.'3 uncer- 
tainty (90% confidence) in the location of the X-ray source. An 
object with an intensity of i<" s =18.4±0.3 mag lies 0''5 from the 
X-ray source. This source is near the detection limit, and chance 
that it is associated with the pulsar s:1.5%. The upper limit on 
the intensity of any object within the error circle is i<" s >18.5. 



ferences between the assumed and measured phases us- 
ing first-and second-order polynomials. The reference 
epochs of the pulse maxima for the two observations 
were 53512.860265(4) and 53539.67325(2) (MJD, TBD). 
The best-fit periods were 10.6112(4) s and 10.6107(1) 
s for 2005 May and June, respectively. This placed a 
limit on the difference in period between the two ob- 
servations of AP < 0.5 ms, or on the period derivative 
of P < 2 x 10~ 10 s s _1 . From the quadratic fits to the 
phases during the individual May and June observations, 
we found AP < 2 and < 0.4 ms (P < 9 x 10" 8 s s _1 
and < 1 x 10~ 8 s s _1 ), respectively. Unfortunately, there 
was a ~2 cycle ambiguity when trying to predict the 
phases over the month interval between observations, so 
the period cannot be refined further using the current ob- 
servations. Nonetheless, the stability of the signal from 
CXO J164710. 2-455216 suggests that it is produced by 
the rotation of neutron star. 

We extracted spectra of CXO J164710.2-455216 using 
standard tools and the acis_extract routine. 13 We found 
no evidence that the intensity of the source changed be- 
tween May and June, so we combined the spectra ob- 
tained from both observations. We mode l ed th e spec- 
trum using XSPEC version 12 (|Arnaudl IT996). All 
uncertainties were computed using la confidence in- 
tervals (Ax 2 = l). The spectrum could be described 
equally well (<10% chance probability) by black body, 
bremsstrahlung, or power law continuum models that 
were absorbed and scattered by the interstellar medium 
(ISM) . To facilitate comparisons with other isolated neu- 
tron stars, we report the parameters of black body model 
(x 2 l v = 63.8/49), for which we derive an absorption col- 
umn of Ah = 1.2 ± 0.1 x 10 22 cm -2 , a best-fit tempera- 
ture of kT = 0.61 ± 0.02 keV, and an apparent radius of 
R hh = 0.27 ± 0.03(L>/5 kpc) km. 

Comparing the three continuum models, we found 
that the interstellar extinction toward the source is 

13 http:/ /www. astro. psu.edu/xray/docs/TARA/ 



A H = (1.2 - 2.6) x 10 22 cm" 2 . This is consistent with 
the range measured for other cluster members, Ah = 
[1.4 — 2.9] x 10 22 cm -2 , where the dispersion can be 
accounted for by variable extinction caused by a fore- 
ground molecular cloud (see Clark et al. 2005). The ob- 
served flux was Fx = 2.41^ x 10~ 13 erg cm" 2 s _1 (0.5- 
8.0 keV). We de-reddened the three models, and found 
L x = 3±2° x 10 33 (L>/(5 kpc) 2 erg s" 1 (0.5-8 keV). 

To evaluate whether CXO J164710.2-455216 is an ac- 
creting system, we searched for an infrared counterpart 
in Aj-band images taken with the Sofl instrument on the 
ESO New Technology Telescope (Fig.EJ). There was no 
counterpart to the pulsar within the 0'.'3 uncertainty in 
its location (90% confidence). The upper limit to its in- 
tensi ty was A s >18 . 5, rul ing out a companion with M>1 
M Q <|Girardi et allEffll . For comparison, the faintest 
stars that have contracted onto the main sequence in 
Westerlund 1 have Mi=2M (Brandner et al., in prep.). 

3. DISCUSSION 

The spectrum and luminosity of 
CXO J164710. 2-455216 (Fig. demonstrate that 
it is not a conventional radio pulsar. First, the power 
lost as the pulsar spins down is too small to produce the 
observed X-ray emission. Assuming that it is a 10 km, 
1.4Af neutron star (/ ~ 10 45 g cm 2 ), the upper limit to 
the period derivative (P < 2 x 10~ 10 s s _1 ) implies that 
the spin-down energy is E = AttIP/P 3 < 3 x 10 33 erg 
s . This is similar to the observed X-ray luminosity, 
whereas magnetospheric emission from a radio pulsar 
with E < 10 35 erg s" 1 would produce L x < 10~ 3 .E 
l|Cheng. Taam. fc Wand l200l . Second, the X-ray 
emission is inconsistent with thermal emission from a 
young, cooling neutron star. Although the characteristic 
temperature of the emission (&T«0.6 keV) is consistent 
with the high en d of the range expected for cooling 
neutron stars fe.g.. lYakovlev fc Pethi ck 2004), the lumi- 
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nosity is 100 times lower than would be expected if the 
surface of the neutron star had a uniform temperature. 

In contrast, the X-ray emission is consistent with 
that of highly-magnetized (B>3 x 10 14 G), slowly ro- 
tating pulsars, which are referred to as magnetars 
l)Duncan fc Thomsponl IT992fl . Known magnetars have 
spin periods between 5 and 12 s, Lx = 1 33 — 10 36 ergs" 1 , 
and spectra that peak at ~2 keV (e.g.. lMereghetti et al.l 
2001 iWoods fc Thompsonll200H) . 

The X-ray emission is also consistent with that from 
faint, wind- accreting neutron star s, such as A 0535+26 
(P=104 s; lOrlandini etaD l200l, X Per (P =837 s; 
iDi Salvo etaLlll998t iDelgado-Marti et alJl2001|) and 4U 
0115+63 (P=3.6 s: ICampana et alJ l200lT However, 
forming a close, accreting binary from a system that 
initially contained two stars with M>40M Q (see be- 
low) and M<1Mq would present a significant challenge 
to models for producing low-mass X-ray binaries (e.g., 
iPorteeies-Zwart. Verbunt. fc Ergmal Il997|) . Therefore, 
we suspect that CXO J164710. 2-455216 is a magnetar. 

It is likely that CXO J164710.2-455216 is associated 
with Westerlund 1, because Galactic X-ray pulsars with 
P=3-30 s are rare (only w20 known; iLiu et al.ll2000|) . To 
evaluate the chance that this pulsar is a random Galactic 
object unassociated with Westerlund 1, we attempted to 
determine the surface density of slow X-ray pulsars on 
the sky. Previous searches for X-ray pulsars were biased 
toward the most luminous X-ray sources, so we searched 
«350 fields in the Galactic plane (— 5<6<5°) that were 
observed with Chandra and XMM-Newton. Aside from 
«15 known examples that were the targets of the rele- 
vant observations, we found no new pulsars with P=3-30 
s (A. Nechita et al., in prep). Therefore, the chance of 
discovering a slow X-ray pulsar in any pointing is <0.3%. 
(We note that the ~15 known slow X-ray pulsars within 
the more conservative latitude range of |6|<1° have a 
surface density of -0.03 degree" 2 . The ACIS-S FOV 
is «0.07 degree" 2 , which also implies that the chance of 
finding a slow pulsar in a Chandra observation is —0.2%.) 
Moreover, CXO J164710. 2-455216 lies only 1.'7 from the 
center of Westerlund 1 (2.3[D/5 kpc] pc in projection). 
Based on the spatial distribution of X-ray sources de- 
scribed above, a source at that location has a —10% 
chance of being a random association. We conclude that 
this pulsar is a member of Westerlund 1 with —99.97% 
confidence. 

The fact that CXO J164710.2-455216 is member of 
Westerlund 1 places a lower limit on the initial mass of 



its progenitor. iClark et"aTI l)2005[) have established that 
the cluster is only 4+1 Myr old, so that the minimum 
mass of a star that could have undergone a supernova 
is «40 Mq. This is supported by the identification of 
several 07-08V stars in the cluster, which have zero-age 
main sequence masses of 34-37 Mq (J. S. Clark et al., in 
prep.). Therefore, CXO J164710.2-455216 was produced 
by a star with M t >40M Q . 

Previously, only three secure lower limits were ob- 
tained on the masses of progenitors to neutron stars, 
all of which were for magnetars. First, a shell of 
HI around IE 1048.1-5937 was interpreted as ISM 
displa ced by the wind of a progenitor with M^=30- 
40M^ llGaensler et alJ 12001 . Second, SCR 1900+14 
was sugge sted to be a me mber of a star cluster <10 
Myr old ( Vrb a et al.l HOOP), placing a lower limit on 
the progenitor mass of M l >20M Q . Finally, SCR 1806- 
20 was claimed to be a member of a s tar cluster that 
is only <4.5 Myr old (although see | Cameron et alJ 
120051 iMcClure-Grifiiths fe Gaenslel I2005D, providing a 
limit of M>50M^fe.g . iFuchs et alJ 11991 Figer et al] 
12005ft . Our evidence that CXO J164710.2-455216 also 
descended from a star with Mi>40 Mq dispels much of 
the doubt that the previous results represented chance 
associations. 

These results demonstrate that some stars with 
Mi>40 Mq do not collapse into black holes at the ends of 
their lives, but instead form neutron stars. This implies 
that massive stars can lose >95% of their mass either 
before or during supernovae. Before super novae, stars 
could lose mass through strong winds (e.g.. iHeger et all 
20031) or be stripped of m ass by binary companions (e.g. 



Wclls tem fe Langeilll999|) . During supernovae, rapidly- 



rotating cores could drive mass away through magneto - 
hydrodynamic winds (e.g., lAkivama fc Wheeled 12005) . 
To determine the importance of these effects, it is neces- 
sary to identify additional neutron stars in star clusters 
and constrain the masses of their progenitors. 
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